Levitation of quantum Hall critical states in a lattice model with spatially correlated disorder 
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The fate of the current carrying states of a quantum Hall system is considered in the situation when the 
disorder strength is increased and the transition from the quantum Hall liquid to the Hall insulator takes place. 
We investigate a two-dimensional lattice model with spatially correlated disorder potentials and calculate the 
density of states and the localization length either by using a recursive Green function method or by direct 
diagonalization in connection with the procedure of level statistics. From the knowledge of the energy and 
disorder dependence of the localization length and the density of states (DOS) of the corresponding Landau 
bands, the movement of the current carrying states in the disorder-energy and disorder-filling-factor plane can 
be traced by tuning the disorder strength. 

We show results for all sub-bands, particularly the traces of the Chern and anti-Chern states as well as the peak 
positions of the DOS. For small disorder strength W we recover the well known weak levitation of the critical 
states, but we also reveal, for larger W, the strong levitation of these states across the Landau gaps without 
merging. We find the behavior to be similar for exponentially, Gaussian, and Lorentzian correlated disorder 
potentials. Our study resolves the discrepancies of previously published work in demonstrating the conflicting 
results to be only special cases of a general lattice model with spatially correlated disorder potentials. 

To test whether the mixing between consecutive Landau bands is the origin of the observed floating, we 
truncate the Hilbert space of our model Hamiltonian and calculate the behavior of the current carrying states 
under these restricted conditions. 



PACS numbers: 72.15.Rn, 71.30.+h 
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I. INTRODUCTION 

A quantum Hall liquid system is believed to undergo a tran- 
sition to a Hall insulator when the strength of the magnetic 
flux density B is turned to zero«i2iiii^ This behavior was sug- 
gested by Khmelnitskii 6 and Laughlir. 7 and explained in terms 
of a levitation scenario where, with decreasing B, the current 
carrying states float up in energy moving one by one across the 
Fermi level. Correspondingly, each time the Hall conductivity 
a xy drops by an amount e 2 /h, until the last current carrying 
state gets depopulated and a xy — > 0. This notion provided 
the basic ingredient for the construction of a global phase dia- 
gram for the quantum Hall effect. 8 In the quantum Hall liquid 
both the longitudinal conductivity a xx and resistivity p xx get 
unmeasurable small in the limit temperature T — > 0, whereas 
the Hall insulator is characterized by a diverging resistivity 
p xx — ► oo while a xx — ► 0. 

However, some of the experimental results are not in accor- 
dance with the floating up picture. For example, direct tran- 
sitions to the Hall insulator from higher (v > 2) quantum 
Hall plateaus have been reported2ii2iii as well as a saturation 
of the Landau level shift as B — > for high quality p-GaAs 
s ample Spi 2 *^ Recently, Yasin and co-workersil have claimed 
that all the conflicting data published so far for different ma- 
terials can be reconciled by a modified global phase diagram. 
They also find that the current carrying states continously float 
up in energy as B — > 0, independent of the appearance of an 
apparent 'metal' -insulator transition at B = 0. 

On the theoretical side, the situation is still unclear. In 
particular, the microscopic origin of the floating remains ob- 
scure because the possible mechanisms put forward hitherto 
are only able to account for a very weak energetical shift, 
but not for the floating of the current carrying states across 
the Landau gap,i&3&i2ii& Therefore, numerical investigations 



were carried out in order to get some hint for the underly- 
ing physical process that causes the critical states to move up- 
ward in energy. Instead of studying the fate of the current 
carrying states when the magnetic field is decreased one can 
also look at the influence of an increasing disorder while B 
is fixed. 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 In both cases, the Landau bands 
eventually start to overlap and the matrix elements connecting 
disorder broadened neighboring Landau bands become impor- 
tant. 

The results of the numerical work for a disordered two- 
dimensional lattice model can be summarized as follows: 
There is a genuine floating of the critical states to higher 
energiesiSLSii^Si 2 ^ ^! This floating is, however, not easy 
to observe in systems with uncorrected random disorder po- 
tentials due to a peculiar annihilation mechanism inherent in 
the lattice model. In this model the Chern states, which cor- 
respond to the critical electronic states in each Landau band 
that are responsible for the integer quantized Hall conductiv- 
ity, get neutralized by the so called anti-Chern states origi- 
nating from the center of the tight-binding band. They move 
outwards in energy almost up to the band edges when the dis- 
order strength is increased. Using spatially correlated random 
disorder potentials instead, it could be shown that the effect of 
the anti-Chern states is reduced so that the expected levitation 
scenario and the floating of current carrying states across the 
Landau gap is seen also in the lattice modeli 2 *^ These results 
also imply that, contrary to other claims, 23 - 26 the direct tran- 
sitions to the Hall insulator occurring within the lattice model 
with uncorrected disorder potentials can of course not be held 
responsible for the direct transitions seemingly observed in 
experiments. 11 Here, finite size effects and limited resolution 
due to finite temperatures are the most probable causes 3 ^ that 
account for the reported behavior. 

The major aims of our work are fourfold: z) We have to 
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check the universality of our previous results which were ob- 
tained for a special model of exponentially correlated disor- 
der potentials. Therefore, we study in the present work the 
influence of Gaussian and, to a lesser extent, Lorentzian like 
correlated disorder potentials in some detail. In particular, the 
results of the Gaussian disorder model will be useful when 
comparing with available analytical 1718 and recent numerical 
work where a 'universal quantitative relation' for the levita- 
tion of extended states has been proposed^i^ 3 .. Also, the cor- 
related Gaussian model may likely become of greater impor- 
tance to future analytical work on this subject because of the 
calculational advantage of Gaussian distributions compared to 
the exponential case. 

z'z) In our previous publication 29 the energetical shift of the 
current carrying states was obtained as a function of disorder 
strength and potential correlation length. However, the ob- 
served shift comprises the broadening effect of the total tight 
binding band and the floating up of the critical states in energy 
relative to the position of the density of states peak. To dis- 
tinguish between these individual contributions and to address 
the question which of the possible definitions for the energet- 
ical shiftASi 2 ^^ are both reasonable and practicable, we have 
to calculate the density of states too. In addition we determine 
the levitation of the critical states as a function of filling factor 
that can be compared with experiments more easily. 

Hi) The general functional dependence of the levitation of 
the critical energy on disorder strength and magnetic field is 
unknown. Besides the conjectures independently suggested 
by Khmelnitski: 6 and Laughlir. 7 there exits only some approx- 
imate analytical work for the weak levitation regimai^iiii^. 
However, even within this limited range the various results 
partially contradict each other. Hence, for the advancement 
of this unsolved problem it seems adjuvant to numerically es- 
tablish some empirical relations between disorder strength W, 
correlation length 77, broadening of Landau levels r and the 
floating up in energy of the critical states. These relations, ob- 
tained for a well defined model, may then serve as a starting 
point for contriving much simpler models and later on as a 
benchmark for checking the outcome of subsequent analytical 
work. 

z'v) To provide a clue for a necessary microscopic expla- 
nation of the levitation scenario, we truncate the disordered 
lattice model and investigate how the floating up of critical 
states changes in the reduced model containing only one or 
two Landau bands. This approach simultaneously provides a 
new method for the scrutiny of the relevant matrix elements. 
This work is already in progress. 



II. MODEL AND METHOD 

The model consists of non-interacting electrons moving in 
the z-y-plane in the presence of disorder and a strong perpen- 
dicular magnetic field. This situation can be described by a 
Hamiltonian defined on the sites of a 2d lattice of size M X L 
with lattice constant a. Choosing the vector potential in the 
Landau gauge, A = (0, — Bx, Q) with (V x A) z = B, we 



have 

(H^)(x,y) = 

w(x, y) ip(x, y) + V [ip(x + a,y)+ ip(x - a, y) 
+ exp(— i2itOiBx/a) ip(x, y + a) 
+ exp(i2naBx/a) ip(x, y — a)] , (1) 

where the magnetic field is chosen to be commensurate with 
the lattice size and as — a 2 eB/h denotes the strength of B 
expressed as the number of flux quanta h/e per plaquette a 2 . 
We mostly use as = 1/8 in the present work and set the units 
of energy and length to V = 1 and a = 1, respectively. The 
on-site disorder potentials w(x,y) are either uncorrected or 
spatially correlated random numbers with probability density 
distributions as described in the next section. 

The density of states p(E, W, rj) has been calculated by di- 
rect diagonalization of square systems of linear size M = 64 
applying periodic boundary conditions in both directions. For 
the calculation of the localization length we utilize a recur- 
sive Green function meihod^i 3 ^ which allows to determine 
the transmission through very long (L ~ 10 7 ) quasi-one- 
dimensional systems of width M in the range 48 < M < 192 
from which the energetical position of the extended states can 
be extracted using a finite size scaling procedure. 

The inverse localization length as a function of energy, dis- 
order strength, and correlation length is defined as the expo- 
nential decay of the trace over the modulus of a sub-matrix of 
the one-electron Green function 

X^(E,W, V ) = - lim ±ln(Tr|G+ L | 2 ). (2) 

L — >oc ALj 

G\ l is the MxM sub-matrix of the Green function G(E) ± = 
(E — H ± ie) _1 acting in the subspace of the columns 1, L 
on the lattice. For the calculation of Xjj(E, W, 7/) periodic 
boundary conditions are applied only in the y-direction (width 
of the system). 

In the case of the projected disordered Harper model we 
used square systems of linear size 48 < M < 96 and ap- 
plied the method of level statistics to determine the energeti- 
cal position of the current carrying states. This information is 
contained in the probability density distribution of the energy 
spacings of neighboring eigenvalues. Previously, this power- 
ful method has been successfully applied to the study of criti- 
cal properties at certain metal-insulator transitions. 36 ' 37,38 ' 39 ' 40 
We have checked that the results calculated by this method 
agree with the ones obtained by the recursive Green function 
method. However, with comparable amount of numerical ef- 
fort, the energetical resolution achieved is better in the Green 
function technique. 

III. CORRELATED DISORDER MODEL 

For the generation of the random disorder potential an algo- 
rithm is required which computes spatially correlated random 
numbers with a known correlation function and an arbitrary 
correlation length, continuously adjustable from zero (the un- 
correlated case) to finite values. To be used efficiently in 
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conjunction with the recursive Green function method which 
operates on narrow but very long stripes, the random poten- 
tial has to be computed "on the fly", i.e., without keeping the 
whole random number field in memory. We generated a corre- 
lated random potential on the lattice Q by averaging over un- 
corrected random numbers e(m) associated with each lattice 
point using a suitable weighting function /. The new corre- 
lated random numbers are defined as 



f(m) 



Yj /(m - m') e(m') 



(3) 



m'ee 



with m = (m X) m y ) <E Q. The disorder potential is w(m) = 
W£(m), where W is the disorder strength and the uncorre- 
lated random numbers e(m) are uniformly distributed over 
[—1,1]. The exact normalization factor TV = X^meO /( m ) 
allows a comparison of different disorder models and cor- 
relation lengths via the homogeneous distribution function 
^(m)( z ) ( see below) and the second moment of the random 
potential. For long-range correlations (77 > 1), N approaches 
the continuum limit (mf for Gaussian correlations), but devi- 
ates considerably for small 77, e.g., r\ < 0.3. 

It is expedient to choose / isotrop, /(m) = /(|m|) and 
vanishing outside some finite circle |m| < R. This cutoff R 
makes the generation of the £ local as we only need to keep 
a window corresponding to an area 2R times the lattice width 
in memory. The £(m) have a local distribution of realizations 
z 6 [—1,1] that, with increasing 77, transforms from a box 
shape into a Gaussian. The form of the distribution function 
can well be approximated by 



P?(m)( Z ) 



N 



erf 



f(0)N- 



4/(0) I Ky/UM-fWN- 



erf ( z 



(4) 



with a second moment < £(m) 2 >,= 172 [/]/3 and g n [f] 
X J2mez 2 /"( m )- The shape of this distribution depends 
on the extent of the weighting function /, hence on the corre- 
lation length, narrowing with longer correlation. This means 
we need a larger W for stronger correlation in order to get 
comparable effects of disorder in the Hamiltonian. 

Choosing /(m) = <d(R — |m|) exp(— |m| 2 /?7 2 ) with suf- 
ficiently large R, we obtain correlated random numbers with 
an approximately (due to the underlying lattice) Gaussian cor- 
relation function 
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exp 



(5) 



with d = m m' and [s] 1 G [0, 1) denoting the fractional 
part of the real number s. The 7(77) are 0.0077, 0.70, and 0.97 
for 77 = 0.3, 0.7, and 1.0, respectively, approaching 7 = 1 
for 77 > 2. For the g 2 [f] we get 1.0, 0.42, 0.17, and 0.04 for 
77 = 0.3, 0.7, 1.0, and 2.0, respectively. 

A Lorentzian like correlated random potential can be con- 
structed by choosing /(m) = 0(i? — |m|) rj 2 /{rj 2 + |m| 2 ). 
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FIG. 1: The peak positions of the density of states p(E, W, rj) in 
the disorder (W) energy (E) plane. The DOS peaks move outwards 
(widening of the total TB band) with increasing disorder for a Gaus- 
sian correlation length 77 = 0.3 (o) and 77 = 1.0 (O). The solid lines 
are quadratic fits SE„ oc — W 2 . 



The corresponding correlation function is 

1 



Jf,(m,m') 



m— m'| 1 

2r} , 



o(|m — m'| 



1 



In this case we have 172 [/] = 2.9 • 10 5 for 77 = 1.0. 



IV. THE BEHAVIOR OF THE DENSITY OF STATES 

Independent of the disorder model a non-zero disorder po- 
tential has two main effects on the density of states (DOS). 
First, the narrow Harper bands that replace the totally de- 
generated Landau levels of the continuum model on the lat- 
tice broaden nearly symmetrically. The total width of a sub- 
band consists of the intrinsic Harper broadening and the ad- 
ditional contribution due to the disorder potentials. For large 
correlation length the disorder induced sub-band broadening 
decreases with increasing Landau level index until the finite 
intrinsic width of the central Harper bands dominates. For 
small correlation length the disorder broadening is nearly in- 
dependent of the Landau level index. The remaining differ- 
ence in the observed sub-band broadening is associated with 
the varying intrinsic width of the unperturbed Harper bands. 
This usually leads to a slightly increasing width of the dis- 
order broadened Harper bands with increasing proximity to 
the tight-binding (TB) band center. The shape of these sub- 
bands depends on the correlation length, becoming Gaussian 
for large 77. Second, with increasing disorder strength the cen- 
ter position of the sub-bands move in energy towards the edges 
of the TB band. This effect is not small and stronger for the 
outer sub-bands. 

The traces of the center positions of the disorder broadened 
Harper bands in the disorder-energy plane are plotted in Fig. [J 
Because of the symmetry of the TB band, only the lower half 
is shown. The widening of the total TB band is proportional to 
the square of the disorder strength within the resolvable range. 
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FIG. 2: The width F of the disorder broadened lowest Harper band 
vs. disorder strength W for different Gaussian correlation lengths 
■q = 0.3 (O), 0.7 (O), 1.0 (A), and 2.0 (O). The solid lines are linear 
fits r oc W. 



The solid lines in the picture are fits W oc y/—SE n for each 
sub-band n, where SE n is the shift in energy of the respective 
DOS peak position due to the disorder W. We note that 6E n 
gets gradually smaller with increasing Landau band index n. 
In contrast, the disorder broadening of each individual Harper 
band is only linear in the disorder strength, r oc W. This 
holds with good accuracy within the accessible range of disor- 
ders, independent of the correlation length as shown in Fig. [2] 
The broadening of the Harper bands is stronger for shorter cor- 
relation length of the disorder potential, however, taking the 
narrowing of the random potential distribution function with 
increasing correlation into account, we find stronger broaden- 
ing for larger correlation length if the second moments of the 
random potential distribution functions are set equal. 

The independence of the observed linearity from the corre- 
lation length enables us to define an effective disorder strength 
W e(W,rj) which causes the same bandwidth E{\V, Tj) — 
r(W e s ) of the lowest Harper band independent of the disor- 
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FIG. 4: The motion of the extended states with increasing disorder in 
the disorder-energy plane is shown for the lowest four Harper bands. 
The states starting at E = -3.3, E = -2, E = -1 are the Chem 
states, while the anti-Chern states start at E ~ —0.2. The Gaussian 
correlation lengths are rj — 0.3 and 0.7, respectively. The dashed 
lines are guides to the eye. 



der model. This measure for the disorder strength allows the 
direct comparison between different correlations length and is 
also readily accessible from experimental data. For a given 
Weff the spreading of the TB band decreases with increas- 
ing correlation length. For large correlation length (?y ^> 4) 
the density of states can be described by the local distribution 
of the random potential folded with the unperturbed Harper 
band structure. For Lorentzian and exponential correlation, 
not shown here, we find qualitatively the the same behavior of 
the DOS. Our results are also in agreement with investigations 
on continuum modelsilii^ 

Due to the linear broadening /\ oc W of the lowest sub- 
band and the quadratic shift 5E\ oc — W 2 of the peak position 
of its DOS as a function of disorder strength, the energeti- 
cal shift is consequently quadratic in the sub-band's width, 
SEi = — /3i(r?) FiiW) 2 . As shown in Fig. [3] the pre-factor 
Pi (rj) behaves for large correlation length like /3i ~ 1 /rj, but 
saturates at a finite value for small correlation length. 



FIG. 3: The inverse of the pre-factor /3i of the quadratic shift, SEi = 
—Pi(rj)r 1 (W) 2 , of the DOS peak position of the first sub-band is 
shown as a function of the correlation length for Gaussian correlated 
disorder and fixed magnetic field B = 1/8. The dashed lines are 
partial linear fits. 



V. THE BEHAVIOR OF THE CRITICAL STATES 

The consideration of finite spatial correlations within the 
random disorder potentials resolves the discrepancy between 
the expectations of levitation of current carrying states within 
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FIG. 5: The same as in Fig[4]with the exception that the correlation 
lengths of the Gaussian correlated disorder potentials are 1.0 and 2.0, 
respectively. 



the continuum model 6 - 7 and the downward moving anti- 
Chern states observed in lattice models with white noise 
disorder42i22i2£iS We will show that both scenarios are only 
two special cases of the lattice model with correlated disorder 
that are associated with different spatial correlation lengths of 
the disorder potential. 

In Figs. |4] and [5] the traces of the critical states for Gaussian 
correlated disorder are shown for several correlation lengths 
?/ in the disorder-energy plane. Again, the upper half of the 
tight-binding band is omitted as it is equal to the lower half 
mirrored at the line E = 0. For short correlation length 
7/ = 0.3, we find no levitation of the lower sub-bands which, 
in fact, move absolutely down in energy following the center 
position of the DOS peaks of the spreading TB band, though 
to a lesser extend. The anti-Chern states which reside near 
the band center for the Harper model without a disorder po- 
tential, rapidly move down in energy with increasing disorder 
strength, annihilating the Chern states from higher to lower 
Landau level index until eventually the last extended state 
vanishes near the band edge. For a slightly larger correla- 
tion length 7/ = 0.7, the lower sub-bands stop moving down- 
wards in energy and remain at their unperturbed energies up 
to a quite large disorder strength. The anti-Chern states still 
move quickly down with increasing W, but less rapidly than 
for smaller ?/. Only shortly before getting annihilated by the 
anti-Chern the extended states from the lower sub-bands start 
to move to higher energy. Now the last extended state is de- 
stroyed slightly above the energy of the lowest unperturbed 
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FIG. 6: The traces of the extended states (Chern and anti-Chern) vs. 
energy for different Gaussian correlation length, r\ = 0.3 (A), 0.7 
(+), 1.0 (o), 2.0 (x), and 4.0 (V) are compared using the effective 
disorder W c s defined by an equal broadening of the lowest Harper 
band. 



Harper band. 

For longer correlation length (see Fig. [3} the extended states 
from the lower sub-bands strongly move up in energy while 
the anti-Chern states still move towards the band edges. The 
energy where the last extended state is destroyed when the 
lowest Chern and the anti-Chern states annihilate moves fur- 
ther towards the band center, passing the unperturbed 2nd 
Harper band energy for rj m 2. If the correlation length is 
increased further the anti-Chern states remain too close to the 
band center to be resolved anymore, while the critical states of 
the lower sub-bands float up in energy clearly across the Lan- 
dau gaps. Thus, for large correlation length, i.e., for smooth 
disorder potentials, we reach the levitation picture anticipated 
from the continuum model. From our data we see no indica- 
tion for a merging of the critical states from the lower sub- 
bands in contrast to another report. 43 

For sufficiently large correlation lengths and disorder 
strengths that produce equal widths of the lowest sub-band, 
the levitation of the 1 st and 2nd critical states becomes slower 
with increasing r\. This behavior is shown in Fig. [6] Also, the 
levitation gets stronger with increasing Landau level index. In 
contrast, the behavior is different for very small 7/ and is prob- 
ably dominated by lattice effects. Similar results have been 
obtained for other disorder models with exponentially 2 ^ and 
Lorentzian like decaying correlation function. 

In Fig. the traces of the extended states' positions with 
respect to the filling factor are plotted. They clearly differ 
from the energy dependence (Figs. |4] and |5Jl since they con- 
tain the behavior of the extended states and the total DOS 
combined. As a function of the filling factor the observed 
levitation is even stronger because of the spreading of the 
tight-binding band and the additional rather small contribu- 
tion originating from the accumulation of the density of states 
at a particular band that comes from sub-bands above (density 
floating 19 - 20 - 24 ). With respect to v, the extended states from 
the lower band always levitate to higher fillings, 31 even in the 
uncorrelated case. The prominent linear slope of the lower ex- 
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FIG. 7: The three figures show the behavior of the extended states 
in the disorder - filling factor plane for Gaussian correlation with 
r\ — 0.3, 0.7, and 1.0. The dashed lines are guides to the eye. 



tended states' traces for not too strong disorder is essentially 
an effect of the DOS behavior. The widening of the TB band 
oc W 2 and the broadening of the Harper sub-bands ex W to- 
gether with either the shift down in energy oc W 2 (rj < 0.3) 
or the absence of any movement (?y > 0.7) leads for a quali- 
tatively Gaussian like shape of the sub-bands to a prominent 
linear term ex W in the filling factor dependence for small W. 
This feature gets less pronounced for larger r\. 

For the disorder model with long-range Lorentzian like cor- 
relations, we find qualitatively the same behavior as for the 
Gaussian correlated random potential. The results are also 
compatible with those obtained previously for exponentially 
correlated disorder^ In Figs. [8] and [9] the traces of the ex- 
tended states are shown vs. energy E and filling factor v, 
respectively, for a correlation length r\ = 1.0. The levitation 
of the extended states from the lowest band across one and a 
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FIG. 8: The behavior of the current carrying states for a long- 
range Lorentzian like correlated background potential with correla- 
tion length rj = 1.0 is shown in the disorder-energy plane. The 
downward movement of the anti-Chern cannot be resolved anymore. 
The dashed lines are guides to the eye. 



half Landau gap in energy and even across two Landau gaps 
w.r.t. filling factor is clearly seen. Also, the second band's 
extended states can be traced nearly across one Landau gap in 
this case. The anti-Chern states from the TB band center are 
moving down too slowly with increasing disorder strength W 
to be resolved in this case. 



VI. DISCUSSION OF THE RESULTS 

The physical problem investigated here is governed by two 
fundamental scales, the energy scale hu c /V, which for the 
lattice model corresponds to Anas, and the magnetic length 
Is/a = (27raB) 1 ^ 2 , where as = p/q is the number of flux 
quanta per lattice plaquette a 2 . For strong magnetic flux den- 
sities B = as{h/ e)a~ 2 , Ib gets smaller than the lattice con- 
stant a, a domain not considered in the present work. Here we 
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FIG. 9: The behavior of the current carrying states for a long- 
range Lorentzian like correlated background potential with correla- 
tion length 77 = 1.0 is shown in the disorder-filling factor plane. The 
dashed lines indicate the linear slope at small W. 
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have chosen 0.3 = 1/8 which guarantees, in connection with 
spatially correlated disorder potentials and our presently avail- 
able computer capacity, the ability to numerically resolve the 
entire ranges of energy and disorder values as shown above. 

To discuss the influence of the correlation length of the dis- 
order potentials one needs a practicable measure that allows 
the comparing of the various results. One choice is the sec- 
ond moment of the disorder potentials which is independent 
of B and can be viewed essentially as fixed for a given physi- 
cal sample. A second possibility is the independent disorder 
broadening of the sub-bands, r oc W, which would be more 
advantageous in comparing with experiments. This method 
has been used in Fig.[6]where traces of the critical states have 
been plotted versus an effective disorder strength. However, 
for weak magnetic field and/or strong disorder this quantity 
may not easily be accessible. 

In considering the observed energetical movement of the 
current carrying states E c , one has to distinguish two contri- 
butions. One is related to the outward shift of the energetic 
position of the individual sub-bands of the DOS with increas- 
ing disorder strength, SE n oc W 2 , which leads to the linear 
dependence on the filling factor v. The second contribution 
is the absolute energetical shift, SE C — E c — E C (W = 0), 
of the critical states (Chern and anti-Chern states) across 
the Landau gaps. Both contributions depend on the correla- 
tion length of the disorder potentials and on the strength of 
the magnetic field. To compare with the energy shift pro- 
posed within the levitation scenario in the continuum limit, 6 - 7 
SE oc (n + 1/2)(w c t)~ 2 , the absolute energetical shift 5E C 
has to be considered. We find a qualitative agreement with 
this relation for larger correlation lengths rj. For comparison 
with experiments, the calculated filling factor dependence of 
SE C (see Figs.0and[5} seems to be better suited. 

The levitation of the energy of the lowest current carrying 
states with respect to the sub-band peak of the density of states 
has recently been proposed to follow the relation^ 

fruj c r\ \Hu! c J 

The definition of an effective shift 5E c s, given by the energy 
difference between position of the DOS peaks and the respec- 
tive current carrying states^i^i is only practical as long as the 
DOS peaks can clearly be resolved which, however, gets in- 
creasingly problematic for small B and strong W. Therefore, 
it is clear, that this definition can, as a matter of principle, 
not be applied for the strong levitation across the Landau gaps 
seen in our investigation. In contrast, the resolution of the 
energetical position of the critical states can continously be 
enhanced by applying larger system sizes. Our results also 
show that even in the weak levitation regime where the en- 
ergetical shift of the critical states SE e g is essentially due to 
the shift of the DOS peak 5E\, the proposed linear relation 
(SEcs)^ 1 ~ rj only holds over a limited range of r\. It strongly 
deviates for small r\ (see Fig. [3} where 8E\ = —Pi(rj)r 2 sat- 
urates. 

In the large disorder regime, where the strong levitation is 
dominated by the absolute energetical shift, even the relation 
5E C r~j W 2 no longer holds. Instead, we clearly observe a 




Energy E 

FIG. 10: The energetical positions of the DOS peaks (□) and the 
positions of the critical states (•) of the lowest two sub-bands for 
the projected disordered Harper model with correlation length 77 = 1 
and magnetic field B — 1/8 are shown vs. energy E. The solid lines 
are the quadratic fits SE oc W 2 . 

dependence with an exponent larger than 2. This behavior 
may change with decreasing B, but, due to the limited system 
size, it was not possible to resolve this problem yet. Finally, 
we did not see any merging of the levitating states, therefore 
true direct transitions are not possible. However, for smaller 
magnetic fields the critical states may get very close so that 
their resolution is not possible any more. 

VII. PROJECTED DISORDERED HARPER MODEL 

In order to get an idea of a possible mechanism that leads to 
the levitation of the critical states, we furthermore investigated 
a disordered Harper model projected onto the lowest few Lan- 
dau levels. This was done by a projection of the highly non- 
diagonal spatial disorder potential within the base spanned by 
the unperturbed Harper eigenstates. Starting by using only the 
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FIG. 11: For the projected disordered Harper model with correlation 
length rj — 1 and magnetic field B — 1/8 the positions of the DOS 
peaks (□) and the positions of the critical states (•) of the lowest two 
sub-bands are shown vs. filling factor v. The solid line is a linear fit 

Sv oc W. 
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FIG. 12: For the projected disordered Harper model with correlation 
length T) = 1 and magnetic field 73 = 1/32 the positions of the DOS 
peaks (□) and the positions of the critical states (•) of the lowest two 
sub-bands are shown vs. energy E. The solid lines are quadratic fits 

SE^W 2 . 
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FIG. 13: For the projected disordered Harper model with correlation 
length rj — 1 and magnetic field B — 1/8 the positions of the DOS 
peaks (□) and the positions of the critical states (•) of the lowest 
three sub-bands are shown vs. energy E. The solid lines are the 
quadratic fits 8E oc W 2 . 



two lowest Harper bands, we could reproduce the widening of 
the total band as well as the levitation of the critical states of 
the lowest sub-band which is shown in Fig. ^|for a Gaus- 
sian correlation length r) = 1. Since there are only two levels 
repelling each other, the center of the 2nd DOS peak moves 
to higher energies, in contrast to the full model. This widen- 
ing leads to a shift of both sub-bands oc W 2 as in the full 
model. The levitation of the lowest critical state is clearly vis- 
ible while the extended states of the 2nd level are essentially 
following the DOS peak. Different from the full model, the 
levitation of the 1st level is with good accuracy oc W 2 over 
the entire resolvable range. The levitation w.r.t. filling fac- 
tor is again linear in W (Fig. II 1>. Further, we do not see any 
floating down of extended states for smaller correlation length 
as it has been observed for instance for Gaussian correlation 
with i] = 0.3 in the full model (see Fig. [4}. 

Dropping the matrix elements describing the inter-band in- 
teraction of the disorder part in the Hamiltonian containing 
only two sub-bands, completely destroys the levitation as well 
as the widening of the total tight-binding band, while preserv- 
ing the localization properties. 

For a smaller magnetic field B = 1/32 the behavior of the 
disordered Harper model projected onto the lowest two Harper 
bands shown in Fig. ^] is qualitatively the same. However, 
the ratio of upward movement of the first sub-band's extended 
states and the downward movement of the corresponding DOS 
peak position changes with magnetic field. The shape of the 
traces remains purely quadratic in any case. 

The addition of the third unperturbed Harper band to the 
projected model adds a qualitatively new behavior which is 
shown in Fig.[0] The DOS peak position of the second band 
remains nearly at the same energy, those of the first and third 
bands are moving outwards, while the extended states from 
the first and second bands move to higher energies and filling 
factors with increasing disorder strength. In the third band, the 
energy of the extended states essentially follows the DOS peak 
position. There is no indication for any absolute levitation of 
these extended states. This corresponds to the behavior of the 



second band in the two-band model shown above. The levi- 
tation of the lowest band's extended states is slowed down by 
the addition of a third band compared to the two-band model. 
Now it resembles closer the trace known from the full disor- 
dered Harper model (Fig. |4}. We conclude, that in order to 
see levitation of a particular Harper band's extended states the 
existence of at least one adjacent band at higher energy is re- 
quired. 

VIII. CONCLUSIONS 

The influence of spatially correlated disorder potentials on 
the levitation of current carrying states has been investigated 
for a quantum Hall system described by a two-dimensional 
lattice model. We considered what happens when the strength 
W of the disorder potentials is increased so that neighboring 
Landau bands start to overlap. Here, we showed our find- 
ings obtained mainly for spatially correlated disorder poten- 
tials which were generated with Gaussian correlations charac- 
terized by a correlation length 77, but results for a Lorentzian 
like correlated disorder have also been presented. They both 
qualitatively agree with the previously obtained outcome for 
exponentially correlated disorden^&S Thus, the results ob- 
tained previously are not the consequence of the special ex- 
ponentially correlated disorder model, but fit into a picture of 
a more general disordered lattice model which for long-range 
correlated disorder potentials agrees qualitatively with what is 
expected from the continuum model without a periodic back- 
ground potential. 

We numerically calculated the density of states p(E, W, 77) 
and studied the scaling of the localization length Am (E, W, 77) 
for various energies E, disorder strengths W, and system 
widths M from which the position of the critical states were 
extracted. Depending on the correlation length, we found 
three different cases. For small correlation length 77 < 0.3 
no levitation of the Chern states, but their annihilation by 
downward moving anti-Chern states has been observed. For 
77 > 2.0 the Chern states float up in energy, at least across one 



9 



Landau gap. No effect of the anti-Chern states could be de- 
tected anymore for 77 > 4 where the anti-Chern remain close 
to the TB band center. In the latter case a levitation scenario 
as expected from the continuum model without a periodic po- 
tential has been observed in our lattice model. In the interme- 
diate case 0.7 < r\ < 1.0 both the levitation of the Chern as 
well as the outwards moving anti-Chern states can be noticed. 
As a function of the filling factor, the shift of the lower crit- 
ical states is always positive, independent of the correlation 
length. One has to bear in mind that the levitation is always 
accompanied by the spreading of the tight-binding (TB) band, 
which is not negligible in particular for small 77. We also have 
to mention that in contrast to a recent report 43 no merging 
of current carrying states was observed in our investigations. 
Therefore, no genuine direct transitions are possible within 
the present model. Of course, also the suggested annihilation 
mechanism of the current carrying states by anti-Chern states 
has to be discarded because this scenario is a special feature 
of the lattice model that cannot serve as an explanation for the 
experimental observations of a seemingly direct transition. 

Besides these results that confirm and complement our pre- 
vious findings for a model with exponentially distributed spa- 
tially correlated disorder potentials, we obtained the following 
new conclusions. i) The introduction of an effective disor- 
der WcfF, which is defined by an equal broadening r of the 
sub-bands, allows a general discussion of the influence on the 
correlation length 77. This is advantageous when comparing 
with other disorder models and experiment. The definition 
was possible only because of the observed linear broadening 
/\ oc W for different Gaussian correlation lengths 77. it) The 
observed movement in energy of the current carrying states 
E c consists of at least two contributions, the shift of the den- 
sity of states peaks 8E n oc —W 2 , and the absolute energetical 
shift which for low disorders follows a relation SE C oc W 2 . 
In the strong levitation regime at larger W, the absolute shift 
SE C and therefore also the observed E c move faster than W 2 . 
The peak of the lowest sub-band follows 6E% oc ~(3\(rf)W 2 
with /3i (77) ~ 1/(3 for rj not too small, but saturates in the 
limit 77 — > 0. Hi) The definition of an effective energetical 
shift of the current carrying states, i.e., 6E e g ~ \Ei — E c \ 
for the lowest sub-band, is practicable only for small disorder 
strength where the DOS peaks can be resolved. Therefore, it 



does neither work in the medium and nor in the strong levi- 
tation regime across the Landau gap as observed in our work. 
Also, the proposed relation (Eq. |6} does not hold for small 77 
where (/3(t7)) _1 is no longer linear in 77. iv) Our study of the 
truncated disordered Harper model, which did not include the 
anti-Chern states, showed the levitation still to exist, provided 
the reduced model contains at least the coupling between the 
sub-band considered and the adjacent sub-band at higher en- 
ergy. The spreading of the TB band as well as the levitation of 
the lowest critical states could already be detected in a model 
projected only onto the two lowest unperturbed Harper bands. 
Dropping the inter-sub-band interaction part of the correlated 
disorder potentials destroyed the levitation and the spreading 
of the TB band while preserving the localization properties. 
These results clearly show that the microscopic origin of the 
critical states' floating up in energy originates in the matrix el- 
ements of the disorder potentials constructed from eigenstates 
of consecutive Landau levels. In the case of the truncated lat- 
tice model the position of the critical states was extracted us- 
ing the method of level statistics. 

To summarize, we presented a comprehensive picture of 
the levitation of the current carrying states in the presence 
of spatially correlated disorder potentials. The behavior can- 
not be described by a simple general relation of the energeti- 
cal shift as a function of magnetic field and disorder strength. 
However, empirical relations have been found which hold in 
certain ranges of the disorder and correlation length param- 
eters. Further investigations, particularly at lower magnetic 
flux densities, are needed to clarify the complete functional 
dependences. Our results for the truncated disordered Haper 
model provide first clues for a microscopical understanding 
of the levitation of the current carrying states at the quantum 
Hall to insulator transition. 
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